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CN . Abstract 

In this paper we derive LPS's criterion for the breakdown of classical solu- 
tions to the incompressible nematic liquid crystal flow, a simplified version of 
Ericksen-Leslie system modeling the hydrodynamic evolution of nematic liquid 
crystals in M 3 . We show that if < T < +oo is the maximal time interval for the 
unique smooth solution u G C°°([0, T), M 3 ), then \u\ + \Vd\ $ L q ([0, T], L P (R 3 )), 
where p and q safisfy the Ladyzhenskaya-Prodi-Serrin's condition: | + 1 = 1 and 
p £ (3, +00] . 
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1 Introduction 

We consider the following hydrodynamic system modeling the flow of liquid crystal 
materials in dimension three (see [2,3,10,12] and references therein): 

u t + u- V«-i/Am + VP = -Ad- Vd, (1.1a) 

d t d + u- Vd = Ad+ |Vd| 2 d, (1.1b) 
V -u = \d\ = 1, (1.1c) 
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for (t,x) G [0, +00) x R 3 . Here u : R 3 — > R 3 represents the velocity field of the 
incompressible viscous fluid, v > is the kinematic viscosity, P : R 3 — > R represents 
the pressure function, and d : R 3 — >■ S 2 represents the macroscopic average of the 
nematic liquid crystal orientation field. We are interested in the Cauchy problem (1.1) 
with the initial value 

(u(0,x),d(0,x)) = (u (x),d (x)) (1.2) 
satisfying the following compatibility condition: 

V -1x0(20=0, \d (x)\ = l, lim d (x) = a e § 2 , (1.3) 

\x\— >oo 

where a is a given unit vector. 

The above system is a simplified version of the Ericksen-Leslie model, which reduces 
to the Ossen-Frank model in the static case, for the hydrodynamics of nematic liquid 
crystals developed during the period of 1958 through 1968 [2,3,10]. It is a macroscopic 
continuum description of the time evolution of the materials under the influence of 
both the flow field u(x,t), the macroscopic description of the microscopic orientation 
configurations d(x,t) of rod-like liquid crystals. Roughly speaking, the system (1.1) is 
a coupling between the non- homogeneous Navier-Stokes equation and the transported 
flow harmonic maps. Due to the physical importance and mathematical challenges, the 
study on nematic liquid crystals has attracted many physicists and mathematicians. 
The mathematical analysis of the liquid crystal flows was initiated by Lin [11], Lin and 
Liu in [12,13]. For any bounded smooth domain in IR 2 , Lin , Lin and Wang [14] have 
proved the global existence of Leray-Hopf type weak solutions to system (1.1) which 
are smooth everywhere except on finitely many time slices (see [5] for the whole space). 
The uniqueness of weak solutions in two dimension was studied by [15,20]. Recently, 
Hong and Xin [6] studied the global existence for general Ericksen-Leslie system in 
dimension two. However, the global existence of weak solutions to the incompressible 
nematic liquid crystal flow equation (1.1) in three dimension with large initial data is 
still an outstanding open question. 

In this paper, we are interested in an optimal characterization on the maximal in- 
terval T that is scaling invariant. So let us first introduce the following definition: 

Definition 1.1. For 1 < p, q < 00, we say a function / = f(t,x) : [0,T] x R 3 — > R is 
in L«([0,T],L p (R 3 )), if 

rp 1 

||/IU«([0,T],LP(M3)) = ( J Q \\f(t, OIIL-fRS)^) 5 , 1 < Q < OO 

= ess sup \\f(t,-)\\ LP{R s), q = oo 
te[o,T] 

is finite. lip = q, then we simply write ||/|| L p([o,t],r») for \\f\\ L p([o,T],LP(R^)- 

We will consider the short time classical solution to (1.1) and address the Ladyzhenskaya- 
Prodi-Serrin's criterion that characterizes the first time finite singular time. The 
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local well-posedness of the Cauchy problem of system (1.1) is rather standard (see 
[5,7,14]). More precisely, if the initial velocity Uq G f/ s (IR 3 ,IR 3 ) with V • Uq = and 
d - a G H S+1 (R 2 ,S 2 ) such that system (1.1) has a unique, classical solution (u, d) in 
[0, T ) x R 3 satisfying 



for any < T < Tq. At present, there is no global-in-time existence theory for classical 
solutions to system (1.1). Thus if we assume T* > is the maximum value such that 
(1.4) holds with T = T*, we would like to characterize such a T*. Motivated by the 
famous work [1], Huang and Wang [7] have obtained a BKM type blow-up criterion 
(see also [16]). However, the techniques involved in this paper are much different from 
[7], which we believe that the result may have its own interest. 

When d is a constant vector field, the system (1.1) becomes an incompressible 
Navier-Stokes equation. Recall that the scaling invariant space L q ([0, T], L P (R 3 )), with 
(p, q) satisfying 

- + - = 1, (1-5) 
p q 

which has played an important role in the regularity issue of Navier-Stokes equation. 
Leray [9] first established the existence of a global weak solution for Navier-Stokes 
equation, now called Leray-Hopf weak solution, that satisfies an energy inequality: 



Although the regularity issue for Leray-Hopf weak solutions of Navier-Stokes equation 
remains open, it is well-known that both uniqueness and smoothness for the class of 
weak solutions of Navier-Stokes equation, in which u G L q ([0, T], L P (R 3 )) for some 
p G (3, +oo] and q G [2,+oo) satisfying Ladyzhenskaya-Prodi-Serrin's condition (1.5) 
have been established through works by Prodi [18], Serrin [19], and Ladyzhenskaya 
[8] in 1960s. On the other hand, for the end point case p = 3, q = +oo , only until 
very recently Escauriaza et al. [4] have finally proved the smoothness for weak solution 
u G L°°([0, T], L 3 (R 3 )) of incompressible Navier-Stokes equation, < T < +oo. 

Motivated by these results for the Navier-Stokes equation, we are going to use 
scaling considerations for system (1.1) to guess which spaces may be critical. We 
observe that system (1.1) is invariant by the following transformation: 



Thus L g ([0, T], L P (R 3 )) is a critical space for (it, Vd) if (p, q) satisfies the Ladyzhenskaya- 
Prodi-Serrin's condition (1.5). 

Our main results are formulated as the following theorem: 

Theorem 1.1. For u G H S (R 3 ) with V • u = and d - a G H S+1 (R 3 ) with 



u G C([0, T), H S (R 3 )) n C\[0, T), H S ~ 2 (R 3 )) and 

d - a G C([0, T), H S+1 (R 3 , § 2 )) n C\[0, T), H S ~ 1 (R 3 , § 2 )) 



(1.4) 




u = lu(l 2 t,lx), P = l 2 P(l 2 t,lx), d = d(l 2 t,lx). 
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I do I = 1 for s > 3. Suppose that (u, d) is a smooth solution to the system (1.1)-(1.2), 
then for given T > 0, (u, d) is smooth up to time T provided that 

IMU<2([o,T],LP(R3)) + ||Vd||L9([0,T],LP(R3)) < + 00, (1.7) 

where (p, q) satisfies the Ladyzhenskaya-Prodi-Serrin's condition (1.5) andp G (3, +00]. 

Notations. We denote by L p , W m,p the usual Lebesgue and Sobolev spaces on IR 3 
and H m = W m ' 2 , with norms || • \\lp, \\ ■ \\w m <p and || • \\h™ respectively. For the sake of 
conciseness, we do not distinguish functional space when scalar-valued or vector- valued 
functions are involved. We denote V = 8 X = (81,82,83), where 8i = d Xi , Vj = 8i and 
put 8 l x f = V 1 f = V(V' -1 /). We assume C be a positive generic constant throughout 
this paper that may vary at different places and the integration domain M 3 will be al- 
ways omitted without any ambiguity. Finally, (•, •) denotes the inner-product in L 2 (M 3 ). 

Remark 1.1. It is standard that the condition (1.3) is preserved by the flow. In 
fact, first notice that the divergence free of the velocity field u can be justified by the 
initial assumption that V • Uq = 0. Indeed, this can be easily and formally observed 
by take V- to the momentum equation. Moreover, applying the maximum principle to 
the equation for \d\ 2 , one also can easily see that \d\ = 1 under the initial assumption 
that I do I = 1- 



2 Proof of Theorem 1.1 

We prove our theorem in this section. Without loss of generality, we assume v — 1. 
The first bright idea to reduce many complicated computations lies in that we just 
need to do the lowest order and highest order energy estimates for the solutions. This 
is motivated by the following observation: 

n/ii 2 ^<^iiav fc /)ii! 2 , v/g^. (2.i) 

This inequality (2.1) can be easily proved by combing Young's inequality and Gagliardo- 
Nirenberg's inequality. 

||V7lk<C(p)||/||£jV fe /HL- \/ feH k (2.2) 
where - - | = -a + (- - |)(1 - a) with i < k. 

p 3 q v r 3 /v ' — 

Now we are in a position to prove our Theorem 1.1. 

Proof of Theorem 1.1 First of all, we note that if p = +oo, Theorem 1.1 has 
been proved in [16], thus let us concentrate on p e (3, +oo). Now for classical solutions 
to (1.1)-(1.2), one has the following basic energy law: 
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Let's concentrate on the case s = 3. For each multi-index a with |a| < 3, by 
applying d£ to (1.1a) and d x +1 to (1.1b), multiplying them by d^u, d x +1 d respectively 
and then integrating them over M 3 , we have 



||||^(ti,vd)lli a + ll^(v«,Ad)|li a 
= -(d%(u ■ v«), <9» - (d%(Ad ■ vd),dy) 



id: 



iQ+l i 



vd),d« +1 d) + (a^ +1 (|vrf| 2 rf),a^ +1 d) 



(2.4) 



i=l 



where I\ a \ ti are the corresponding terms in the above equation which will be estimated 
as follows. Now for \a\ — 1 in (2.4), integrating by parts and using the divergence free 
condition V • u — and (2.2), we arrive at 



|/ lfl | = \(dl(u-Wu),dlu)\ = mu-Vu,d l x u)\ < Cj R3 \Vu\*dx 

3p 18 2p 



<c\\ u \\i + /\\v 2 u\\ L t p <c\\ u \\r P 3 \\u\\r P 3 



Mv 2 u\\h- 



Combining Cauchy's inequality, Sobolev's inequality and the fact |V<i| 2 
(since \d\ = 1) gives 



(2.5) 
-d ■ Ad 



= \(d x (Ad-Vd),d x u)\ = \(Ad-Vd, Au) | 

< C\\Ad\\l 3 + -L\\V 2 u\ 

3p 18 

<c\\vd\\it"\\^d\\ii p 

< C\\Vd\\p\\Vd\\t 3 + U\W 3 d\\li + HV 2 W ||| 2 ). 



|2 

<C\\Vd\\W\\V 3 d\\^ + ±\\V 2 u\\h 



(2.6) 



Similarly, 



= | (d 2 x u ■ Vd + d x u ■ Vd x d, d 2 x d) \ 

< C\\Ad\\^ + l 2 \\V 2 u\\l 2 + C\\Vuf L3 

2p p 2p p 

I p— 3 1 1 1 1 p— 3 I II Y~7 J|| p— 3 || v~7 J|| p— 3 1 



< C'(|| M ||£; 3 |h||^ 3 + ||Vrf||£; 3 ||Vrf||i P , 
+U\\V 3 d\\l 2 + \\V 2 u\\h), 

= \(dl(\Vd\ 2 d),d x Ad)\ 

<C(||Vrf|| 6 L6 + ||V 2 d||! 3 ) + ^||VAd||| 2 

< C\\V 2 d\\l + l 2 \\VAd\\ 2 L2 

^cuvdnl^iivrfii^ + ^iiv^iii,. 



(2.7) 



(2.8) 



Taking the above estimates (2.5)-(2.8) in (2.4) for |a| = 1, we arrive at 



|||(V^,V 2 d)||| 2 + C||(VVV 3 d)||| 2 

< c(\\ u \\i- 3 \\u\\t 3 + \m\i~ 3 \\vd\\t 3 ). 



(2.9) 



Next we derive an estimate for \\u\\ p LP and ||Vd||^ p . First of all, we multiply (1.1a) 
by \u\ p ~ 2 u and integrate over M 3 to obtain that 



pdtW a \\LP 



j w ,\ur 2 \Vu\ 2 dx+\{p-2)j R 
= -(VP, \u\ p - 2 u) - (AdVd, \u\ p - 2 u). 



\u 



p— 4 



\V(\u\ 2 )\ 2 dx 



(2.10) 



Observe that 

Ad- Vd = V- (Vd© Vd- ^|Vd| 2 /), (2.11) 

where Vd © Vd denotes the 3x3 matrix whose the entry is given by <9jd • <9jd 

for 1 < i,j < 3. And taking div to (1.1a), we arrive at 

AP = -divdiv(u © « + Vd Vd - ]-\Vd\ 2 I). 

An application of the //-estimate of elliptic systems to the above equation, there exists 
P(t) such that 

J R3 \P-P{t)\^dx 

< Cj R3 \u\ p+2 + \Vd\ p+2 dx 

< C{(J R3 \u\vdx) E T L (J^\u\^dx^ 
+(f R3 \Vd\vdx) E ^(f RS \Vd\^dx) l v} 

< C{(J R3 \u\ p dx)^(J R3 \u\ p - 2 \Vu\ 2 dx)i 
+(/ R8 |Vd|'dx)^(/ R8 |Vdr 2 |V 2 d| 2 d*)f} 

< c{(/ K3 H%)^ + (/ K3 |vw^} 
+ULsJur 2 \vu\ 2 dx +j R3 |vdp- 2 |v 2 d| 2 dx) 

< G{ii^#hii p L P + iivdi#iivdir LP } 

\ur 2 \Vu\ 2 dx + / R , |VdMV 2 d| 2 dx). 

Thus we have 

| (VP, \u\ p ~ 2 u)\ = \(P - P(t),div(\u\ p - 2 u)}\ 
< C J R3 \V\u\ 2 \\u\ p ~ 4 \u\\P - P(t)\dx 



< j R3 ^ivi«i 2 i 2 | M r 4 + c\u\p+ 2 + \p- p{t)\^dx 

< c{hi#ii* + iivdi#iivdir LP } 

+^(/r3 |VM 2 |>r 4 + \u\ p ~ 2 \Vu\ 2 dx + J R3 \Vd\P- 2 \V 2 d\ 2 dx). 
Similarly, applying (2.11) and by Sobolev's inequality and Cauchy's inequality, we get 
| (Ad- Vd, \u\ p ~ 2 u)\ = \ (div(VdQVd-l\Vd\ 2 I),\u\ p - 2 u}\ 

< c-diuii^iHiip + nvdii^iivdiiip} 

|V|u|TM p ~ 4 + \u\ p ~ 2 \Vu\ 2 dx + / R3 |Vdp- 2 |V 2 d| 2 dx). 
Putting the above two inequalities into (2.10) we have 



dt\\ U \\LP 

2p 



+ / R8 \u\ p - 2 \Vu\ 2 dx + \ J R3 \u\ p -*\V(\u\ 2 )\ 2 dx 



_2e_ ZJK ' ' ' Vl ' " (2.12) 

< c{||HI£;1l^ll! P + l|vd||£-||vd|r LP } + ^/ R3 |vdr 2 |v 2 d| 2 dx. 

Next differentiating (1.1b) with respect to x, we have 

d x dt - Ad x d = d x (\Vd\ 2 d - u ■ Vd). (2.13) 
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We multiply (2.13) by |Vd| p 2 d x d and integrate over M 3 to obtain that 

\i\\Vd\\% + /„ \Vdr^d\ 2 dx+ \{p- 2)/ R3 |VrfMa,(|Vd| 2 )| 2 rfx 
= (|Vd| p - 2 x d,0 x (|Vd| 2 d-ii- Vd)) 

= || Vd|| p + 2 2 + (| Wl'-^d, 2WdWd x dd) + (V • (| Vdr 2 ^d), « • Vd> 

< C(|| Vdjj£& + ||u||g? a ) + La |Vdr 2 |V 2 d| 2 d* < 2 - 14 ) 

< C-dltill^lHI^ + UVdll^llVdll^} 
+s(Jr" |«MV«| 2 ds + /„ |Vd|^ 2 |V 2 d| 2 dx). 

Combining (2.12) and (2.14) gives 

+ HVd|| p LP ) + /„ |«r 2 |V«| 2 + |Vdr 2 |V 2 d| 2 da; 

_2£_ (2.15 

< C{||«||i;'|H|^ + ||Vd||£;»||Vd||S,}. 

2 2 

Now by (1.5) we have g = thus ||w||£; 3 and || Vd||£~ 3 belong to L^O, T]. For p, 
we divide into two case: 

Case 1. If 4 < p < +oo, then ^ < p. Thus combining (2.9) and (2.15) gives 

|(ll (v«,v 2 d) ||| 2 + 1| («,vd) ||£ P ) + c\\ (vy v 3 d) ||| 2 

< c(i| M |#i| W ir LP + iivd|#nvdir LP + i). 

By Gronwall's inequality, we have 

||(v M ,v 2 d)(t,.)|li2 + ||( M ,vd)(t,-)|| p iP + c£||(vVv 3 d)( s ,-)|li2d s 

2p 2p 2.16 

< CeMCf^\Hs,.)\\r p 3 + ||Vd( S ,-)||£; 3 + l)ds] < +oo, 
for all t E [0,T]. 

Case 2. If 3 < p < 4. Multiplying (2.15) by (\\u\\ p LP + ||Vd||^ p )^, we obtain 

i(\\u\\ p LP + nvdir LP )^ < c(HH# + iivd#)(iMr LP + iivdf LP )^. (2.17) 

We add (2.17) to (2.9) to obtain 

liii(v^v 2 d)iii 2 + anr LP + uvdii^] + cikw v^iii. 

< c(|H# + ||vd#)(|| w ||i P + ||vd|r LP )^. 

By Gronwall's inequality, we get 

||(v«, v 2 d)(t, .)||| 2 + (\Ht, oil?, + IW, onw^ 1 

+C/J||(VVV 3 d)( S ,0||i2d S <C'exp[C'/ c ;||( M ,Vd)( S ,.)||£; 3 d S ]<+oo, 
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for all t e [0,T]. 

Next for \a\ =3. For 7 3jl , we need to use the following Moser-type inequality (see 
[17, p. 43]): 

\\D s {fg)\\v < C(|MM|V S /IU 2 + ll/IU-IIVffl^). (2.19) 



Thus we have 

I4il 



= \{d%.div(u <S> u), d£u)\ 
<C\\V\u®u)\\i 2 + MV*u\\l 2 

< C\\u\\U\VMh + U^u\\h (2.20) 
<C{||V«||l 2 ||V 4 M ||I 2 ||V M ||l 2 ||V 4 M ||l 2 } 2 + ^||V 4 M ||i 2 

<c\\vu\\% + i\\vM\h- 

For I 3j2 , we apply (2.11) and (2.19) to obtain that 

\h,2\ = mv ■ {Vd®vd-\\vd\ 2 i),dtu)\ 

< C\\V\Vd®Vd-\\Vd\ 2 I)\\\ 2 + \\\V A u\\ 2 L2 
<C||Vrf||| oc ||V 4 rf||| 2 + |||V 4 W ||| 2 (2.21) 

< C{||V 2 d||l 2 ||V 5 c/||| 2 ||V 2 rf||I 2 ||V 5 rf||| 2 } 2 + i\\V*u\\l 2 
<G||V 2 rf||I 4 2 + |(||V^||| 2 + ||V 5 rf|| 2 2 ). 

Similar in the proof of (2.21), I 3 3 , 1 3 i can be bounded as follows: 

tasi =\mu-vd),d*d)\ 
<c\mu-vd)\\i 2 + ±\\v 5 d\\i 2 

< C(||Vd||i.||V s «||i a + IMIUV 4 d||| 2 ) + ^||V 5 d||i 2 

<c{||v 2 d||l 2 ||v 5 d||i 2 ||v M ||| 2 ||vS||| 2 (2 - 22) 

+ l|V W ||I 2 ||V 4 M ||l 2 ||V 2 rf||l 2 ||V 5 rf||l 2 } 2 + ^||V 5 d||i 2 

< c(\\vw£ + liv^llii) + |(||v 4 M ||| 2 + ||v 5 rf||| 2 ), 

and 

\hA =\m\vd\*d),d 5 x d}\ 
<c\\dmvd\id)\\i 2 + i 2 \\vsd\\i* 

< C(||V 3 |Vrf| 2 ||i 2 + |||Vd| 2 ||i„||V s d||£ a ) + 3^11 V 5 d||i 2 

< C(||Vd|| 2 .||V 4 d||i 2 + ||Ad||UV 3 d||i 2 ) + ^| V 5 d||| 2 ( 2 - 23 ) 

< C(HV 2 d||i 4 + {||V 2 d||| 2 ||V 5 d||l 2 ||V 2 d||l 2 ||V 5 d||| 2 } 2 ) + ^||V 5 d||i 2 
<C||V 2 d||i 4 + |||V 5 d||i 2 . 

Putting (2.20)-(2.23) into (2.4) for \a\ = 3 and by (2.16) and (2.18), we arrive at 

|(l|v 3 HI! 2 + llv 4 d|il 2 ) + ||v 4 w ||i 2 + ||v 5 d||i 2 

<C(\\V 2 d\\% + \\Vu\\}*) <+oo. 

Integrating the above inequality with respect to time from to t G [0, T] and by (2.3), 
we conclude u(T, •) e H 3 (R 3 ) and d(T, •) - a e H 4 (R 3 ). Thus the proof of Theorem 
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1.1 is completed. 
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